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ABSTRACT 


The  basic  formulation  of  scattering  problems  in  terms  of 
integral  equations  is  examined,  for  the  special  case  of  perfectly 
conducting  bodies  of  revolution.  In  particular,  the  singularities 
of  the  integrals  which  arise  in  this  context  are  studied,  in 
relation  to  the  transition  to  a  numerical  solution  by  means  of 
the  method  of  moments.  For  both  H-field  and  E-field  equations,  it 
is  found  that  finite  matrix  elements  can  be  deduced  in  a  way  that 
is  uniquely  determined  by  the  integrals  themselves.  No  ad  hoc 
procedures  are  required  to  secure  convergence,  but  one  such 
procedure,  which  is  commonly  used,  is  shown  to  be  capable  of 
giving  accurate  results  for  those  integrals  which  tend  to  diverge 
logarithmically.  In  the  H-field  solution,  other  integrals  also 
arise,  which  yield  finite  terms  not  normally  included  in  the 
theory,  when  the  necessary  limiting  procedures  are  carried  out. 
These  terms  can  play  a  significant  role  in  the  case  of  a  body 
whose  profile  contains  a  corner.  -  _ 
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INTRODUCTION 


The  formulation  of  scattering  problems  in  terms  of 
integral  equations  is  a  well-established  procedure,  and  the 
application  of  numerical  techniques  to  the  solution  of  these 
equations  has  spawned  a  large  and  varied  literature.  The 
integrals  in  these  formulations  are  generally  improper,  in  the 
sense  that  the  integrands,  being  essentially  Green's  functions, 
have  integrable  singularities.  The  integrals  must  therefore  be 
treated  with  some  care,  particularly  when  the  transition  is  made 
to  a  set  of  linear  equations  by  means  of  the  method  of  moments. 
The  diagonal  elements  of  the  resulting  matrices  are  directly 
related  to  regions  of  integration  where  these  singularities 
manifest  themselves.  Many  analyses  have  been  carried  out  in  which 
divergent  integrals  are  replaced  by  finite  expressions  or 
convergent  integrals  in  a  somewhat  ad  hoc  way,  the  procedure 
being  justified  by  physical  reasoning.  It  is  always  recognized 
that  these  divergences  arise  because  of  oversimplifications  in 
the  application  of  the  numerical  procedure,  such  as  the  use  of 
delta-function  expansion  or  test  functions  in  the  method  of 
moments.  It  is  assumed  that  errors  due  to  the  approximations  used 
will  turn  out  to  be  small,  and  indeed  this  is  so  when  the  results 
are  compared  to  known  cases,  such  as  the  sphere.  However,  one 
never  knows,  a  priori,  that  the  calculations  will  be  as  accurate 
when  the  technique  is  applied  to  a  body  of  radically  different 


shape.  It  is  generally  accepted  that  the  method  of  moments 
procedure  itself  will  converge  more  rapidly  (as  a  function  of  the 
number  of  points  used  to  represent  the  body)  if  its  matrix 
elements  are  computed  with  increased  precision,  and  it  may  even 
be  possible  to  gain  assurance  that  the  answers  to  which  the 
method  converges  are  in  fact  the  correct  ones.  It  therefore  seems 
desirable  to  formulate  the  problem  in  a  way  that  represents  the 
integrals  accurately  during  the  transition  to  a  system  of  linear 
equations . 

In  this  note  the  singular  kernels  are  examined  in  the 
special  case  of  scattering  by  conducting  bodies  of  revolution, 
and  it  is  found  that  finite  diagonal  matrix  elements  can  be 
deduced  in  a  way  that  is  uniquely  determined  by  the  integrals 
themselves.  Although  the  expressions  found  are  slightly  different 
than  those  customarily  used  in  the  literature,  it  is  not  expected 
that  the  numerical  results  which  follow  from  their  application 
will  be  very  different,  with  the  possible  exception  of  bodies 
with  corners.  In  this  case  the  present  evaluations  (for  the 
H-field  solution)  contain  terms  not  generally  used  and  which  are 
directly  related  to  the  abrupt  change  of  the  tangent  direction  of 
the  profile  of  the  body  at  a  corner. 

The  analysis  is  carried  out  here  in  terms  of  the  H-field 
equation,  since  it  is  an  integral  equation  of  the  second  kind 
with  no  derivatives  operating  on  the  unknown  function.  For  these 
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reasons,  it  has  the  better  chance  of  numerical  stability,  at 
least  away  from  resonant  frequencies  of  the  corresponding 
internal  problem.  In  the  usual  treatment  of  the  E-field  problem, 
derivatives  of  the  expansion  and  test  functions  appear,  which  are 
approximated  by  difference  operations.  This  adds  another  element 
of  inaccuracy  to  the  analysis,  producing  errors  which  are 
difficult  to  assess  in  advance.  Here  also,  divergent  integrals 
occur,  and  these  can  easily  be  treated  by  the  methods  given 
here.  The  E-field  solution  is  discussed  in  the  final  section  of 
this  report. 
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1.  FORMULATION  of  the  PROBLEM. 

Let  S  be  a  closed  surface,  covering  a  perfectly  conducting 
body,  and  let  H0(r)  be  the  magnetic  field  which  would  be 
produced  by  a  system  of  sources  in  the  absence  of  this  body.  With 
the  body  present,  the  total  field  produced  by  this  same  system  of 
sources  satisfies  the  equation: 


H(r)  =  H  (r)  +  ]  K(r')  *  V'g(r,r')  dS’ 
°  S 


where  r'  is  the  variable  of  integration  on  the  surface,  V'  is  the 
gradient  relative  to  r',  dS '  is  the  surface  element,  K  =  n  *  H, 
g(r,r')  is  the  scalar  Green's  function: 


,  R  -  r-r 


ikR 

g(r,r')  *  -ywjf 


and  k  is  the  wavenumber.  The  suppressed  time  dependence  is  of  the 
form  and  only  unit  vectors,  indicated  by  an  underline, 

are  specifically  identified  as  vectors  in  this  report.  The  field 
point,  r,  is  exterior  to  the  surface  in  the  above  equation  and  K 
depends  on  the  unknown  values  of  tangential  H  on  the  surface.  By 
letting  r  approach  the  surface,  and  taking  the  limiting  value  of 
the  integral,  we  obtain  the  Maue t ^  equation  for  K: 

K  (r)  =  i  K ( r )  -  n  *  J  K(r')  *  V'g(r,r')  dS ' 

U  ^  T  « 


where  K0  =  £  *  H0.  In  this  integrand,  r  is  on  S,  and  the 
surface  integral  is  actually  convergent.  The  formulation  of  the 
scattering  problem  in  terms  of  this  equation  for  a  body  of 
revolution  is  well  known, [2]  but  we  give  the  basic  steps  here 
in  order  to  make  the  development  reasonably  self-contained. 

We  define  the  function 

G  (R )  =  { 1-ikR)  SLi*i 


so  that 

V'g(r,r ' )  =  <r-r' )  G(R ) 
and  the  Maue  equation  becomes 

K  (r)  =  i  K ( r )  +  J  n  *  [(r-r')  *  K(r' ) ]  G(R)  dS'. 
°  1  S 


We  let  £  represent  either  of  two  orthogonal  tangential  unit 
vectors  at  r,  £'  being  one  of  the  unit  vectors  at  r',  and  then 
write  the  integral  equation  in  component  form.  Indicating  vector 
components  by  superscripts,  we  have: 

Kp(r)  -  l  Kp(r)  +  l  J  Mpp'(r,r')  Kp(r’)  dS ' 

°  p'  S 

where  l  sums  over  the  two  values  of  p',  and  the  dyadic  kernel  has 


components 


MPP  ( r , r ' )  =  £  *  n  *  [(r-r‘)  *  £' ]  G ( R )  . 


Next  we  specialize  the  surface  to  that  of  a  body  of 
revolution,  taking  the  axis  of  symmetry  as  our  z-axis,  and 
introduce  a  coordinate  system  for  points  and  vectors  on  S.  The 
three-dimensional  cylindrical  (P,<t>,z)  system  will  also  be 
employed,  where  P  is  the  distance  from  the  z-axis  and  is  the 
azimuthal  angle.  We  suppose  that  the  body  half-profile  is 
defined  by  the  relation 
p  =  p(z) 

and  introduce  the  arc-length,  t,  along  the  body  by  means  of  its 
differential  element: 

dt  =  (dz2  +  dP2)1/2. 


The  derivative,  dP/dz,  is  related  to  8,  the  angle  which  the 
tangent  to  the  profile  makes  with  the  z-axis,  by  the  equation: 


dP  _ 

3?  ' 


tan 


8 


and  in  terms  of  8  we  have 

dz  =  cos  8  dt  ,  dP  =  sin8  dt  . 

For  our  purpose  it  is  useful  to  think  of  t  as  the 


independent  variable,  measured  from  a  reference  point 
(z0»p0)  on  the  body,  and  to  express  8  as  a  function  of  t. 


Then  z  and  P  are  obtained  as  integrals: 


t 

z  =  zQ  +  J  cosB(s)  ds 
t 

p  =  pq  +  j  sin3(s)  ds. 
fco 

These  equations  will  play  a  basic  role  in  the  analysis  of  the 
singular  portions  of  the  integral  equation  for  K.  The  coordinates 
of  a  surface  point  are  now  t  and  ♦  ,  and  the  surface  element 
itself  is  dS  =  Pdtd<j>. 

Let  t  be  a  tangential  unit  vector  in  the  direction  of 
increasing  t,  and  let  ♦  be  a  similar  vector,  in  the  ( P,  -plane , 
in  the  direction  of  increasing  <J>.  Together  with  the  normal  unit 
vector,  n,  these  form  a  local  orthogonal  system.  We  assume  that 
the  postive  t-direction  is  defined  such  that 
n  =  ±  *  t  . 

It  is  useful  to  have  the  notation 
£(£)  =  £  *  £ 

so  that 

£(  t )  =  ,  £(  jM  =  -t.  . 

Returning  to  the  dyadic,  M,  we  have 

Mpp  (r,r‘)  =  [(r-r1)  x  £’  ]  ’  {£  x  n)  G(R) 
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=  £  •  ; 'r-r1 )  X  £'  ]  G (R ) 

=  {(£  x  r )  •£'  +  £  •  (£’  x  r '  )  }  G (R ) 


The  dot  products  can  be  conveniently  evaluated  in  the  cylindrical 
coordinate  system,  and  this  is  carried  out  in  Appendix  1  for  the 
more  general  case  in  which  r  is  off  the  surface  a  small  distance 
in  the  normal  direction.  When  specialized  to  a  surface  point  for 
r,  we  obtain  the  dyadic  elements  in  the  form  derived  by 
Harrington  1 21 s 


Mtt(r,r')  =  { [ ( p'-p)cos S'  -  (z'-z)  sin8']  cos^ 

-  2  p  cos  6'  sin2(V2)}  G (R ) 

Mt<*>(r,r')  =  (z'-z)  sint  G(R) 

M^t(r,r')  =  [p'sin8  cos  8'  -  P  cos  8  sin8'  -  (z’-z)  sin8  sin  8'  ]  • 
sin>P  G (R ) 


M’y(r,r') 


{[(p'-p)cos8  -  (z'-z)  sin8]  cos'!* 
+  2  p'  cos  8  sin2(V2)}  G(R) 


In  these  formulas,  ♦ 

The  dyadic  depends  on  the  azimuthal  angles  only  through 
as  a  consequence  of  the  cylindrical  symmetry,  and  the  next  step 
is  the  modal  expansion.  Substituting 

Kp(r)  =  I  Kp(t;m)  eim * 
m 


8 


and  defining 


KS(t;m)  =  J2  kS  e"lra*d4> 

O  4«  g  O 


and 


LPP'(t,f;m)  -  j2’  MPP'(r,r')  elm*  d* 

0 


we  obtain  the  uncoupled  modal  integral  equations 

Kp(t;m)  =  -j-  Kp(t;m)  +  l  J  Lpp  (t,t';m)  Kp  (t'jmjp'dt'. 

p‘ 


In  R,  the  argument  of  G,  the  introduction  of  cylindrical 
coordinates  yields 


R  =  I  r-r '  =  [(p'-p)2  +  (z'-z)2  +  4  p  '  p  sin2 ( <|>/2  )  ]1/2  . 


Since  R  is  an  even  function  of  ♦  ,  the  ♦-integrals  involved 
in  the  components  of  L  reduce  to  the  following: 


Ti(t,t':m)  =2  J  G(R)  cos^  cos(m4»)  d^ 
1  0 


T_(t,t';m)  =2  J  G(R)  sin2('|»/2)  cos(m<P)  d'P. 
£  0 


IT 

T,(t,t,;m)  =  2  i  )  G (R)  sin'l'  sin(m'l')  dt. 
*  0 
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In  terms  of  these  integrals  we  have  the  desired  expressions  for 
L,  again  in  Harrington's  form: 

Ltfc  =  [(P'-P)  cosB'  -  (z'-z)  sinB']  -  2p  cos  3'  T2 

Lt<(>=  (z'-z)  T3 

L^t=  [p'sinB  cos3'  -  P  cos  3  sin3'  -  (z'-z)  sin  3  sin  3'  ]  T3 
[(p'-p)  cos  3  -  (z'-z)  sin3]  +  2  p'  cos  3  T2  . 

In  the  simplest  application  of  the  method  of  moments,  the 
unknown  functions  KP(tjm)  are  expanded  in  pulse  functions  and 
the  error  resulting  from  their  use  is  forced  to  zero  at  discrete 
points.  The  pulse  functions  are  each  unity  over  short  intervals 
of  t,  and  zero  elsewhere.  They  are  non-overlapping  and  together 
cover  the  whole  range  of  t.  The  test  points  lie  within  these 
intervals,  not  necessarily  at  the  mid-points.  We  are  therefore 
required  to  evaluate  integrals  of  the  type 

^2+^2  , 

]  Lpp' (t. ,t' )  p'  dt' 

fc2-6l 

where  t^  is  an  arbitrary  test  point  and  the  integral 
corresponds  to  a  pulse  function  whose  non-vanishing  range  covers 
another  test  point,  t2»  (From  now  on  the  mode  number  is  not 


retained  explicitly  in  the  notation).  Since  the  quantities  LPP' 
will  be  obtained  by  numerical  integration  over  4>,  it  is 
impossible  to  carry  out  the  subsequent  t-integrals  analytically, 
and  hence  the  integral  above  is  further  approximated  by  the 
simple  expression 

Lpp,'t1,t2)  P2  A 

Here  P2  is  the  value  of  P  at  t2  and  A  =  6^  +  62.  If  one 
wished  to  use  more  general  test  functions,  the  t-integrals  would 
in  any  case  have  to  be  evaluated  numerically,  and  one  would  again 
be  led  to  the  consideration  of  integrals  of  L  alone  over  short 
t-segments,  as  above. 

This  simple  approximation  will  not  suffice  when  ti  * 
t2 ,  because  the  ^-integrals  defining  the  components  of  L(t,t') 
all  fail  to  converge  when  t  and  t'  coincide.  These  singularities 
are  integrable,  however,  and  by  considering  in  detail  the 
behavior  of  L(t,t')  when  t'  *  t  is  small,  we  can  obtain  finite 
expressions  for  the  integrals  desired  in  the  case  t2  =  tj. 


2.  ANALYSIS  of  the  SINGULARITIES. 

In  the  expressions  for  the  components  of  the  dyadic  M, 
given  in  the  previous  section,  many  of  the  terms  vanish  when  t  = 
t',  since  this  implies  that  z  =  z'  and  P  -  P'.  To  determine  the 
behavior  of  these  terms  when  t*  -  t  is  small,  we  return  to  the 
formulas  for  z  and  P  as  integrals  over  t.  It  is  now  assumed  that 
the  profile  of  the  body  is  such  that  the  tangent  angle,  B,  is 
continuous  with  continuous  derivatives,  near  t,  so  that  the 
expansion 

1  .  •  2 

0*  =  B+0T  +  ie  t  +  ... 

is  possible,  where 
t  =  t'  -  t 

and  the  dots  refer  to  derivatives  with  respect  to  t.  The 
profile  may  well  include  corners,  but  they  are  modelled  as  smooth 
curves  with  very  small  radii  of  curvature.  When  this  series  is 
substituted  in  the  integrals  for  z  and  P,  we  obtain  the 
expansions : 

z'  *  z  +  cos 8  T  -  -j-  B  sin  B  t2  +  ... 

1  2 

p'  =  P  +  sinB  t  +  -^  6  cos  B  t  +  ... 
where  all  unprimed  quantities  are  evaluated  at  t.  From  these 


formulas  and  the  expansion  for  3'  itself,  we  obtain 

1*2  3 

(p'-p)  cos  3  -  (z'-z)  sinB  =  j  3  t  +  0(t  ) 

(p'-p)  cos  3’  -  (z'-z)  sin  3’  =  -  ^  3  t2  +  0(t3) 

and  . 

2  2  2  1  2  4 

(p'-p)2  +  (z'-z)2  =  T  "  "T2  (B)  t  +  **• 

Therefore,  when  t'  -  t  is  small,  the  formulas  for  the  components 
of  L  reduce  to  the  following: 

•  • 

Ltt=  (-  3  t2+  ...)  Tl(t)  -  2 p  ( cos 3  -  3  sin 3  t  +  ..)  T2(  t) 

Lt4>=  (cos 3  t  +  ...)  T 3 (  t ) 

L4’t=  (-3  P  t  +  ...  )  T3(  T) 

(|  3  :2  +  ...)  Tl(T)  +  2  ( p  +  sin3  t  +  ...)  cos  3  T2(t)  . 

The  higher  order  terms,  represented  by  the  ellipses  in  these 

expressions,  will  not  contribute  to  the  t '-integrals,  to  the 
order  of  accuracy  of  interest,  as  will  be  seen  below. 

The  T-integrals  depend  on  t  and  t',  and  hence  on  t, 
through  R.  The  t-dependence  has  been  recognized  in  the  notation 
above,  but  the  fixed  parameter,  t,  is  suppressed.  We  make  the 
definitions 
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J-^l — L 


r  -r-  r  -i  j  v  \  t  •  rj'  v~^~- 


"7  »"J  «T.” 


'  V  J"  ■".■ 


X*  =  (p'-p)2  +  (Z'-Z)2  *  T2  +  ... 


a 


=  p'p  =  p  (p+  sinS  t  +  ...) 


and  substitute  in  the  formula  found  earlier  for  R  to  obtain 


2  =  x2  +  4a2  sin2(V2). 


R 


Thus,  for  x  =  0  and  t  very  small,  x  is  zero  and  R  is  proportional 
to  which  is  the  cause  of  the  failure  of  convergence  of  the 
T-integrals  at  this  point. 

To  deal  with  these  integrals,  we  recall  the  definition  of 
G  and  note  that 

(1-z)  ez=  (l-z)(l+2  +  ^  z2  +  ...)  -  1  “  j  z2  ~  \  +  ••• 

We  can  therefore  put 


K  ( l-z)ez  = 


+  H  (  z  ) 


2z 


where 


H(z)  =  j  [(l-z)eZ  -1  +  ^~2  ] 
z  J 


and  express  G  in  the  desired  form: 


-  I  -  1  + 

3  8  * 
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ft*  ■■  «  I 


G(R)  =  {-J— - - 

( ikRj  2ikR 


- —  +  H  (  ikR )  } 


G<R>  -  TT-  Mr  +  -To  +  h(R)  ] 


where  h  is  defined  by 

h(R)  =  ( ik ) 3  H(ikR)  . 

The  new  quantity,  h(R),  is  a  bounded  function  of  R  with  the  value 
ik3 

— 3 —  at  R=0 .  We  substitute  this  representation  of  G  in  the 

T-integrals  to  find 

,  n  2 

T ^  ]  [— ■*"  'Jr  +  h(R)  ]  cos  cos(m^)  d^ 

0  R 

l  it  .  ,2 

t2  =  -J7  i  [-3  +  |r  +  h(R)]  sin2(*/2)  cos (m ♦)  d* 

0  R 


r  1  .  k2 


T3  =  -~j  }  [— =-j  +  -Jr  +  h(R)  ]  simf*  sin(m'l')  d*  . 

0  R 


The  terms  of  these  integrals  involving  h  are  convergent  for  any 
value  of  t,  and  the  singularities  obviously  arise  from  the 
inverse  powers  of  R. 

The  next  step  is  the  expression  of  the  remaining  factors 
in  each  of  the  integrands  as  the  sum  of  two  terms,  one  of  which 
yields  a  convergent  integral  even  when  multiplied  by  R~^,  while 


I 


the  other  allows  the  singular  portions  to  be  evaluated  in  closed 
form.  The  key  to  this  separation  is  the  possibility  of  a  change 
of  the  variable  of  integration,  to  a  new  quantity: 
u  =  2  sin(t/2). 

Since 

du  =  cos (t/2)  dt  , 

we  note  that  integrals  involving  a  power  of  R,  the  factor 
cos(t/2)  and  a  simple  function  of  sin(t/2)  will  be  integrable  in 
closed  form.  Starting  with  T]_#  we  define  the  quantity 

2 

A^t)  =  cost  cos(mt)  -  cos(t/2)  [l  -  sin2(t/2)]  , 


and  write 


cost  cos(mt) 


2 

cos(  t/2)  [l  -  ^±1  sin2(  t/2)  ]  +  A^  t)  . 


The  sin2 ( t/2 )-term  has  been  supplied  and  its  coefficient  so 
chosen  as  to  make  Ai(t)  proportional  to  t^  as  t  tends  to 
zero.  Then  the  integral  of  Aj/R^  is  convergent  even  at  t  = 

0,  and  the  divergent  portions  of  Tj  are  completely  contained  in 
the  integral 


U, 


=  -4^r  1  (~r  +  cos (  t/2)  [l  - 


4m2+ 


i  i  n  2  (  t/2)  ]  dt. 


~>  TT 


The  new  variable  is  now  introduced,  and  becomes 


U1  = 


2" 


2 

i  { 


0 


.  2  2 ,  2,3/2 
(a  u  +x  )  / 


2 (a2u2+  x2)1/2 


Id 


4m2  +  3 


u  )  du . 


The  term  involving  u2/R  is  convergent  when  x=0 ,  and  it  is 
removed  to  be  lumped  with  the  other  benign  integrals.  Since 


2  2  2  2 
u  du  _  1  I  du  x^  |  du 

_  .  2  2  2,3/2  '  2  J  .22,  2.1/2  ~  2  ,  2  2^  2,3/2 

0  (a  u  +x  )  '  a  0  (a  u  +x  )  '  a  0  (a  u  +x  )  '  , 


i 


we  need  only  the  integrals 


i  du  1  ,  r2a  +  (4a2  +x2)^2  i 

J  — 5  2 - 2~T72  =  a  log  L - 1 - m -  J 

0  (a2u2+x2)i/2  a  Xf 


and 


,  du 

nJ  .22^  2,3/2 
0  ( a  u  +x  ) 


2, a  2  2,1/2 

c  (4a  +x  )  ' 


to  complete  the  evaluation  of  U^.  Collecting  these  steps,  we 
have  shown  that 


TT 


i 

0 


h(R)  ]  {cos (  w 2)  (i-im_±l 


sin  2  (  */2  )  ] 


+  A,  ( 'I' )  }  d  'P 


where 


U,  = 


1 

2  tt 


2 

J  ( 


,  4m  +3  2 

1 - s —  u 


2 

k/2 


0 


(a2u2+x2)3/2  (a2u2+x2)1^2 


}  du 


1  ,,  ,  4m2+3  x2,  ,,  2,  2  * -1/2 
— j  (1  +  — g —  — j)  (4a  +  x  ) 

fx  a 


4ira' 


1  ,  4m  +3  .  2  2 ,  , 

__  ( - - k  a  )  log 


2a  +  i 

!  4a  +x  ) 

x  j 

and 

it  2 

=  2  J  {g  (R )  cos 'i'  cos(m'l')  -  cos (  t/2 ) 

2 

- ^-cos(V2)  [l  -  4lV3  sin2(  »/2)  ]}  d». 

4ttrj  * 


The  integral  is  convergent  at  t=0,  and  its  integral  is 
well-behaved  for  all  t.  For  small  values  of  x,  we  have  the 
expansion 


U1  = 


2Tta  x' 


4  Tia 


, 4m2+  3  .22.,  ,4a 

( - - - k  a  )  log  ( 


FI 


)  + 


which  can  be  used  to  assess  the  diagonal  matrix  element  terms 
containing  Tj.  Only  the  leading  term  of  this  expansion  will  be 
needed  in  the  end.  It  should  be  noted  that  P,  and  hence  a,  is 
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assumed  not  to  be  small,  so  that  a  thin  wire  is  excluded  from  our 
analysis . 

The  other  two  integrals  are  treated  in  the  same  way, 
beginning  with  the  definitions 

A2(40  =  sin2(V2)  [cos(m'J')  -  cos  (  ip/2  )  ] 

A3(<f0  =  sin<l>  [sin(m^)  -  2m  sin(if»/2)] 

=  sin  ^  sin  (m  <(')  -  4m  cos  (  <J»/2  )  sin2(V2) 


These  A-functions  are  again  0('l>4),  and  T2  and  T3  are 
written 


T2  =  2 TT"  J  f  ^3  +  +  h(R)  ]{cos(  i>/2)  sin2  (  V2 )  +  A„  (  ♦)  }  d  * 

0  R  1 


n-  i  2 

T3  =  2J~  I  f  ~3  +  20  +  h(R)  3  {4m  cos(  ty/2 )  sin2(  ♦/ 2)  +  A^  (  '(')  }  d  ^ 

0  RJ  4 


The  only  terms  of  these  integrals  which  are  divergent  at  t=o  are 


11  •? 

U2  =  T?  3  — r  cos('l'/2)  sin  ('|>/2)  d ^ 

0 


0J  <a2u2+x2>: 


) 


=  4im  U2  . 


This  integral  has  already  been  evaluated,  and  to  complete  the 


analogy  to  T^,  we  define 


T  2  —  U  2  +  W 2 


T3  =  U3  +  W3  , 


where 


U  =  — log 

8*aJ 


(4a2+  x2)-1/2 

4  "a 


U3  =  4im  U2 


W_  =  2  J  {g  (R )  sin  (  <P/2  )  cos(m'l') - ~  cos  (  i>/2  )  sin  iV2)}  d  Y 

0  4TTR-5 


2i  ]  {G(R)  sin  ^  sin  (m  <l>) 

0 


—^3  cos  (  V2  )  sin  ('ly'2)}  d  'P 

TTR 


Like  Wj ,  W2  and  W3  are  well-behaved  and  convergent  at  t=o, 


3.  EVALUATION  of  the  MATRIX  ELEMENTS. 


As  mentioned  in  Section  1,  to  apply  the  method  of  moments 
we  need  to  evaluate  integrals  of  the  form 

t2+  2  , 

I  LPP  (t.,f)  P'df 

t  -<S  1 

Z2  1 

where  t^  and  t2  are  test  points  of  the  procedure.  When  these 
points  are  distinct,  the  simple  approximation 

Lpp'(t1,t2)  P2  (51+52> 

will  suffice,  and  the  T-integrals  required  are  obtained  by 
numerical  integration.  The  original  integrands  of  these  integrals 
will  usually  be  sufficiently  well  behaved  to  cause  no  difficulty 
here.  The  mimimum  value  of  R  attained  in  such  an  integral  is  just 
x,  and  in  this  case 

X  =  l(P2~P1)2  +  {z  2-z  l)2]l//2 

is  the  chord  connecting  the  profile  points  t^  and  t2»  The 
decomposition  of  the  T-integrals  into  the  explicitly  evaluated 
U-functions  and  the  well  behaved  W-integrals  would  be  a  useful 
aid  to  numerical  integration  if  kx  were  small,  but  this  is  not 
likely  to  be  the  case  unless  an  excessive  number  of  points  had 
been  used  to  describe  the  body. 


Turning  to  the  diagonal  matrix  elements,  we  write  the 


required  integrals  over  t*  in  the  form 
52  , 

J  LPP  (T)  p'dt 

-51 

where  x  =  t'  -  t  as  before,  and  the  t-dependence  is  suppressed  in 
the  integrand.  In  evaluating  these  integrals  we  shall  retain  only 
terms  of  the  order  A  =  6^  +  62  or  larger.  The  formulas  of 
the  preceding  Section  for  LPP '  are  now  restated: 

Ltt  =  -  \  0  t2  (u1+w1)  -  2 P  cosP  (u2+w2) 

Lt<*>  =  cos0  x  (u3+W3) 

*t 

=  -  6  P  t  (U3+W3) 

=  \  0  j2  (U1+W1)  +  2  P  cos  0  (U2+W2)  . 

In  Section  2  the  exact  formulas  for  the  U-functions  were 
expressed  in  terms  of  x,  which  in  the  present  context  is 

x  =  [(P'-P)2  +  (z'-z)2]1/2  =  x  -  -L.  (S)2  t3  +  ... 


It  is  easily  verified  that  x  may  be  replaced  by  t,  and  the 
leading  terms  of  the  expansions  of  these  functions  used  to 
evaluate  the  integrals  over  t  to  the  required  accuracy. 
Similarly,  a  may  be  replaced  by  p  in  these  expressions,  which 
then  become  simply 


U 


1 


U 


2 


■ - -y -  T  •  •  • 

2*ptz 

ITT"  1109  'f¥ 


1}  + 


•  •  • 


and 

U3  =  4im  U2 

The  W-integrals,  of  course,  are  0(1)  with  regard  to  t,  and  hence 
the  contributions  of  Wi  and  W3  to  the  final  matrix  elements 
which  contain  them  are  of  lower  order  than  A.  since 

j  t  log  t  di  «  J  t2  (log  t  -  ^  )  , 

even  U3  fails  to  survive  the  integration  over  t,  and  thus  $ 
and  L't’t  are  zero  to  the  desired  accuracy. 

We  are  left  with 

Ltfc  =  --|bt2ui-2  p  cos  6  (U2+W2) 


2 


and 
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in  which  we  have  again  replaced  the  factor  p'  by  P.  In  addition 
W2<t)  has  been  replaced  by  its  value  at  t  *  0,  since  higher 
terms  in  the  series  expansion  of  W2  would  not  contribute.  We 
express  the  logarithmic  integral  in  the  following  way: 

62 

j  log  ( |  t|  )  dT  =  «2  log(62/e)  +  6^  log(61/e) 

=  uA  log(-H|)  +  (l-w)A  log  [  ] 

where  u  is  defined  by  the  relations 

A 2  - 

=  (l-W)A. 

Then 

tJ  )  d  T 

=  A  log(A/e)  +  A  [y  logy  +  (1-y)  log(l-y)] 


where 


r 


(1-tO 


(1-u) 


The  value  of  the  parameter,  r,  depends  on  the  position  of 
the  test  point  in  question  within  the  range  of  the  pulse  function 
which  contains  it.  In  general,  the  profile  of  the  body  is  defined 
by  a  sequence  of  pairs  (zm,Pm),  and  there  is  no  guarantee 
that  test  points  will  fall  at  the  mid-points  of  pulse  function 
ranges  unless  this  is  deliberately  arranged.  If  u  happens  to  be 
one-half,  then  r  is  also  one-half,  and  in  any  case 

7  -  r  -  1  * 

Collecting  results,  we  have  obtained  the  following  formula 
for  the  diagonal  matrix  element: 


j  Lfct( t)  p»  d  t  =  - 


-6. 


A  8 


+  Wr  4  Uoo(-SA-)  +  i  -  8«p3  w,(o>] 


4e  P 


and,  of  course,  =  -  Ltfc 


In  the  integral  over  +  which  defines  W2(0),  we  have 


simply 


R  =  2p  sin ( ty/2 ) 


Using  the  definitions  of  W2  and  G(R),  we  can  write 


8.P3  W2(0)  -  I  I  lltl  ifnTVJT  ' 


where 

F('l')  =  (1  -  ikR)  exp(ikR)  cos(m'l')  -  cos(V2) 


The  complete  integrand  vanishes  at  <|>  =  0,  and  for  the  purpose  of 
numerical  evaluation,  the  leading  term  of  the  expansion 


F  (  ♦) 

sinTF/ST 


(k2P2  - 


m 


■)  *  + 


may  be  employed  for  small 


4.  DISCUSSION  of  the  H-F I ELD  SOLUTION 


Our  results  for  the  diagonal  matrix  elements  differ 
from  those  of  Harrington because  of  the  treatment  of  the 
divergent  integrals.  Harrington  avoids  divergence  of  all  the 
T-integrals  by  moving  the  test  point  slightly  off  the  surface,  by 
a  distance  equal  to  A/4  in  our  notation.  Diagonal  matrix  elements 
are  then  evaluated  by  the  same  approximation  as  off-diagonals, 
but  the  integrals  are  now  finite  even  for  ti  =  t2»  In  the 
formulas  for  LPP ' ,  the  coefficients  of  all  the  T-integrals 
except  T2  vanish  in  this  case,  yielding  =  0,  but 

the  T^-term  in  Ltfc  and  is  lost,  and  the  T2  contrib¬ 

utions  are  evaluated  somewhat  differently.  We  defer  discussion  of 
the  terms  ±A6/4tt,  which  are  obtained  from  Tj,  and  show  first 
that  the  two  evaluations  of  the  T2~terms  are  in  fact  very 
similar . 

Aside  from  terms  that  vanish  with  x,  we  have  shown  that 
T2(T)  =  "  lo3<|  T|}  +  Ko 

8itpJ  1  1  0 

where  K0  is  a  constant,  evaluated  in  Section  3.  When  we 
integrated  over  x,  we  obtained 
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to  this  accuracy.  If  we  had  placed  the  field  point,  r,  off  the 
surface  by  a  distance  C  along  the  positive  normal,  we  would  have 
found 

R  =  [x2  +  4a2  sin2(  V2)l  1/2 

as  before,  but  where  x  and  a  now  have  the  more  general 
definitions 

x2  =  (P’-P)2  +  (z'-z)2 

-2?[(p'-p)  cosB  -  (z'-z)  sinB]  +  ?2 

and 

a2  =  p'  (  p+C  cos  B)  . 
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< 


I  . 

t  \- 

I 


If  now,  we  let  t  -  t',  we  obtain  simply 

2  2 
x  =  ? 

a2  =  P( P  +  c  cos  3)  , 

and  since  C  is  small,  the  asymptotic  expansion  derived  in  Section 
3  is  valid  for  T2 r  but  now  in  terms  of  x.  Again  neglecting 
terms  which  vanish  with  C,  the  value  of  T2>  with  T  =  0  and  C 
small,  is  the  same  as  the  expression  above,  but  with  C  replacing 
t.  Consequently,  if  C  is  chosen  equal  to  r^/e,  and  the  T2 
integral  evaluated  numerically,  this  procedure  should  give  the 
same  results  as  our  formulas,  to  the  desired  accuracy.  Recalling 
the  range  of  r,  we  see  that 


which  covers  the  value  1/4  used  by  Harrington. 

Regarding  the  new  term,  ±AS/4n,  we  combine  it  with  the 
diagonal  contribution  of  1/2  made  by  the  explicit  term 

j  Kp(t;m) 
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in  the  modal  integral  equation.  This  combination  is  then 


1  /l  _  AfS.  t 

7  (1  K 


in  the  t-equation,  and 


7  d  +  7?>  K* 


in  the  4>-equation.  For  a  convex  portion  of  the  body  profile,  A3 
is  negative  by  our  sign  conventions.  For  a  smooth  surface  and  a 
reasonable  density  of  points,  A3/2*  will  be  small  compared  to 
unity,  and  the  new  term  will  have  little  effect.  It  may,  however, 
be  comparable  to  the  remaining  part  of  Ltfc  or  L^.  The 
comparison  here  is  between  A3  and 

cos 3  [log  (-^)  +  K'  ]  , 


where 


K'  =  1  +  log (-5-)  -  8wp3  W-(0)  . 


The  ratio  A/p  will  generally  be  small  (thin  wires  having  been 
excluded),  and  the  expression 


A  ,  .  A  . 

Ip  1o9  <pM 


I 


can  not  exceed  1/e,  for  values  of  A  less  than  P.  Thus  Ag  should 


be  compared  to  the  bound 
cos 8  (i  +  K'  |) 


for  the  remaining  terms  of  this  matrix  element. 

The  most  likely  situation  in  which  the  new  term  could  have 
an  effect  is  in  the  immediate  vicinity  of  a  corner  of  the 
profile,  as  at  the  base  of  a  cone.  In  order  to  accomodate  this 
term,  the  scheme  for  describing  the  body  must  allow  for  the 
pre-processing  of  a  sequence  of  profile  points  to  yield  pulse 
function  boundaries,  test  points  and  values  of  S  or  Ag  at  the 
latter.  If  a  pulse  function  spans  a  corner,  the  corresponding  Ag 
should  equal  or  closely  approximate  the  actual  change  in  tangent 
direction  at  this  point.  In  this  case,  the  total  diagonal  matrix 
elements  can  be  significantly  altered. 

The  combination  of  the  Ag  terms  with  the  existing  1/2  in 
the  diagonals  yields  expressions  which  are  very  similar  to  the 
formulas  given  by  Poggio  and  Miller for  the  H-field  integral 
equation,  when  the  field  point  falls  at  a  corner.  In  our 
notation,  their  term  is 


1_ 

2 


K? 
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which  applies  to  either  value  of  p.  We  assumed  that  the  surface 
was  actually  smooth,  although  the  radius  of  curvature  at  a 
'corner'  of  the  profile  could  be  extremely  small  compared  to  the 
wavelength.  Since  the  limiting  process  implicit  in  the  derivation 
of  the  Maue  equation  involves  taking  the  field  point  right  up  to 
the  surface,  there  is  no  natural  scale  against  which  non-zero 
radii  of  curvature  may  be  compared.  In  the  case  of  the  formula  of 
Poggio  and  Miller,  a  mathematically  precise  corner  is  assumed, 
and  their  modified  term  applies  only  when  the  field  point  falls 
at  this  corner.  For  the  t-equation,  our  formulation  will  give  the 
same  result  as  theirs,  but  the  formulas  differ  for  the 
^-equation.  This  is  less  disturbing  than  it  might  be,  due  to  the 
fact  that  the  ^-component  of  surface  current  is  (integrably) 
infinite  at  a  mathematically  precise  corner.  Maue m  makes  the 
case  for  this  statement  in  a  general  context,  and  he  obtains  the 
same  behavior  for  the  component  of  current  along  an  edge,  such  as 
the  junction  of  a  cone  with  a  flat  base,  as  is  exhibited  by  the 
exact  solution  for  a  wedge. 

Returning  to  the  first  equation  given  in  Section  1,  we 
know  from  potential  theory  that  the  tangential  component  of  the 
surface  integral,  as  a  function  of  r,  is  discontinuous  at  the 
surface,  suffering  a  jump  equal  to  K(r)  at  this  point.  We  also 
know  that  when  r  is  placed  on  the  surface,  this  integral  is 


convergent,  in  the  sense  that  the  limit  approached,  when  a  small 
region  covering  r  is  exluded  from  the  integral,  is  independent  of 
the  shape  of  this  region,  as  its  maximum  chord  goes  to  zero. 
Maue's  equaton  results  from  the  fact  that  the  value  to  which  this 
integral  converges  is  the  average  of  its  values  when  approached 
from  opposite  sides  of  the  surface.  If  we  did  not  know  this,  we 
could  still  obtain  an  integral  equaton  for  K  by  actually  taking 
the  limit  of  the  integral 

n  x  j  K( r ' )  *  V1 g(r ,r ' )  dS ' 

S 

as  r  approached  S  from  the  outside.  As  an  exercise,  this 
procedure  is  carried  out  in  Appendix  2  for  the  individual  modal 
integral  equations  with  the  result,  of  course,  that  the  term 
KP ( t ;m )  is  replaced  by  1/2  KP(t;m),  and  the  remaining 
integrals  all  converge.  The  validity  of  this  derivation  depends 
entirely  on  the  integral  Ti  and  the  asymptotic  expansion 
derived  for  it  in  Section  3.  The  general  technique  may  be  useful 
in  other  situatons,  such  as  the  form  of  the  E-field  equation 
discussed  in  Section  5,  where  a  similar  surface  integral  is 
involved  which  can  be  evaluated  in  finite  form  as  a  limit  (as  r 
approaches  S),  but  in  which  one  cannot  put  r  actually  on  S  and 
still  have  a  convergent  integral. 
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5.  THE  E-FIELD  SOLUTION 


If  E0  is  a  source  field,  and  r  is  a  point  exterior  to 
the  surface  S  of  a  perfect  conductor,  then  the  electric  field 
satisfies  the  equation 

E(r)  =  E  (r)  +  i  moo  j{K(r')  g(r,r')  +  [K(r')‘V']  V'g(r  ,r  • )  }dS ' 

s  li 

where  K  is  again  n  x  H  on  the  surface  and  w  is  the  magnetic 
permeability.  This  form,  of  course,  follows  directly  from  the  use 
of  the  standard  dyadic  Green's  function^4),  and  the  second  term 
of  the  integrand  is  usually  transformed  as  follows: 

J  [K  ( r  '  )  •  V*  ]  V'g  ( r  ,r  '  )  dS  • 
s 

=  -  v  J  K(r')  •  V'g(r ,r ' )  dS ' 
s 

-  v  J  g  (  r  ,r  '  )  v  •  K  ( r '  )  dS  '  . 

The  field  point  is  now  allowed  to  approach  the  surface, 
and  the  boundary  condition  n  *  E  =  0  is  applied,  with  the  result 

n  x  E  (r)  =  Lim  n  *  j  g(r,r')  K(r')  dS ' 
r  +S  S 

+  hr  Lim  (nxV)  J  g ( r ,r ' )  V'  •  K(r')  dS  ’ 


From  potential  theory  we  know  that  this  surface  integral  and  its 
tangential  derivatives  are  continuous  functions  of  r,  on  or  off 
the  surface,  hence  the  limits  are  taken  by  simply  placing  r  on  S 
and  no  convergence  difficulty  arises. 

The  method  of  moments  is  now  applied ^1,  with  expansion 
and  test  functions  in  the  Galerkin  procedure.  The  derivative  of 
the  expansion  function,  inside  the  integral  over  r',  is 
approximated  by  a  difference  operation,  and  the  V  outside  the 
integral  is  transferred  to  the  test  function  upon  integration 
over  r,  where  it  is  also  approximated  by  a  difference. 

When  carried  out  in  detail,  the  following  'J'-integrals 

occur : 


TT 

Si  =  2  ]  g(R)  cos(mf)  d  ♦ 

1  0 


7T 

S -  =  2  j  g(R)  cos'}*  cos(m'J')  dt 
0 


TT 

S-j  =  2i  ]  g(R)  sin'}'  sin(m'l')  d  'I' 


and 


The  leading  term  in  the  expansion  of  g(R)  for  small  R  is  simply 
1/4  irR,  hence  S3  is  convergent  and  the  others  are  simple.  The 
divergent  portion  of  Si  and  S2  is  just 


L.  s'  a* 

"  A  K 


1  < 2  du 

2"  0  ( a2u2  +  x2)1/2 


This  integral  is  logarithmic  and  it  can  be  dealt  with  just  as 
T2  was  in  Section  3.  In  particular,  the  field  point  can  be 
removed  a  distance  ri/e  from  the  surface  and  all  terms  become 
finite.  No  contribution  involving  A0  arises.  This  is  consistent 
with  the  E-field  equation  of  Poggio  and  Miller,  since  the  term  in 
their  equation  affected  by  the  presence  of  a  corner  vanishes  when 
the  boundary  condition  for  a  conductor  is  applied. 

Another  E-field  equation  can  be  obtained  from  the  first 
equation  of  this  section  by  letting  r  approach  the  surface  and 
then  equating  the  tangential  components  of  E  to  zero: 

— n  *  E  (r)  =  Lim  n  *  J  {K(r)  g(r,r') 

r  +S  S 

+  Kr  [K  ( r  '  )  •  V]  V'g  ( r  ,  r  '  )  }  dS  •  . 


The  differential  operators  are  now  treated  exactly,  since  the 
derivatives  are  carried  out  explicitly  on  the  scalar  Green's 
function,  at  the  cost  of  a  more  complicated  kernel  for  the 
integral  equation  (see  Andreasen ^ 1 ) .  The  integral  will  be 
convergent  only  if  r  approaches  S  after  the  integrations  are 
performed . 

The  methods  of  this  report  have  been  applied  to  this 
equation  and  it  is  found  that  all  integrals  remain  finite,  but 
the  analysis  is  much  more  complicated  in  detail.  Nine  'I'-integrals 
occur,  many  with  inverse  fifth  powers  of  R,  and  the  asymptotic 
forms  are  more  tedious  to  obtain.  A  particular  feature  arises 
which  may  well  cause  difficulty  in  a  numerical  analysis  based  on 
these  equations,  however.  One  of  the  diagonal  matrix  element 
terms  turns  out  to  be  inversely  proportional  to  A,  the  pulse 
function  range  in  the  t-direction.  If  our  surface  integral  is 
used  to  find  the  electric  field  at  a  point  close  to  a  plane 
surface,  due  to  a  small  region  rear  the  point,  an  expression  is 
obtained  which  is  inversely  proportional  to  the  linear  dimensions 
of  that  region,  as  the  point  approaches  the  surface.  The  formula 
obtained  is  strongly  dependent  on  the  shape  of  the  region,  and  in 
our  case  is  proportional  to  1/A.  The  problem  is  not  one  of 
convergence,  since  A  is  not  infinitesimal,  but  it  may  lead  to 
numerical  instability,  thus  defeating  the  purpose  of  this  scheme. 


APPENDIX  1.  EVALUATION  of  the  DYADIC  COMPONENTS 


Let  z,  £_  and  J  be  unit  vectors  in  a  cylindrical  coordinate 
system.  The  directions  of  _P  and  depend,  of  course,  on  location, 
and  for  two  points  on  the  surface  of  the  body,  the  relations 

p  •  p '  =  4>  *  <$>'  =  cos  ^ 


and 

•£'=-£•  _$•  =  s in  ^ 

are  valid,  where  f  =  <P’  -  <t>.  The  normal,  n,  and  the  tangent  unit 
vector,  _t,  are  expressed  in  terms  of  z  and  _P  by  the  equations 

t.  =  cos3  z  +  sin3  _P  ,  n  =  -sin3  z  +  cos3  _P 

If  the  field  point,  r,  is  off  the  surface,  in  the  direction 
of  _n,  by  a  small  amount  C,  then  we  can  write 

r  =  zz  +  P  P  +  Cn=  (z  -  C  sin  3)  z  +  (p  +  t  cos  3)  P  . 


Working  out  the  cross-products  such  as  z  x  £  =  <f>,  we  find  that 
_t  x  r  =  (p  cos  3  -zsinS  +  O 

and 


4>  x  r  =  (z  -  t  sin3)  p  -  (p  +  ?  cos3)  z 


At  r',  the  point  is  on  the  surface,  hence  the  analogous 


APPENDIX  2.  DERIVATION  of  the  MODAL  EQUATIONS  as  LIMITS. 

If  the  analysis  of  the  main  report  is  repeated  with  the 
field  point,  r,  off  the  surface  by  a  small  distance  C  in  the 
direction  of  the  positive  normal,  the  modal  equations  will  be 

Kp(t;m)  =  Kp(t;m)  +  Lim  I  j  Lpp  (t,t';m)  Kp  (t';m)  p'dt' 

C+o  p’ 

Using  the  results  of  Appendix  1,  one  finds  that  the  dyadic 
components  are  now  given  by 

Ltfc  =  [(p'-p)  cosB'  -  (z'-z)  sinB'  -  C  cos(B'-B)]  T^ 

-  2  (p+C  cos  3)  cos  6'  T 2 

Lt<f)  =  (z'-z  +  ?  sin  6)  T3 

L<*>t  =  [  p  *  sinB  cos  B'  -  p  cos  B  sinB'  -  (z'-z)  sin  B  sinB' 

-  ?  sinB']  T3 

and 

L<*m*’  =  [(p'-p)  cosB  -  (z'-z)  sinB  -  C]  T^ 

+  2  P '  cos  B  T0 


The  T-integrals  have  the  same  definitions  as  before,  but 
as  noted  in  Section  4,  we  now  have 

R2  =  x2  +  4a2  sin2  ( ty/2 )  , 

where 

x2  =  (P’-P)2  +  (z'-z)2 

-  2  ?  [  ( p  *  —  p )  cos  6  -  (z'-z)  sinB]  +  C2 

and 

a2  =  p'  ( p+e  cos  8) 

The  entire  discussion  of  these  integrals  is  still  applicable, 
including  the  asymptotic  expansions  for  small  values  of  x. 

We  can  simply  put  C  =  0  everywhere,  except  in  the  vicinity 
of  the  point  t'  =  t,  where  we  must  evaluate  the  effect  of  the 
C-terms.  In  those  terms  of  the  dyadic  components  which  are 
identical  to  terms  already  discussed  (with  r  on  the  surface),  we 
can  immediately  put  C  =  0  for  any  t',  since  these  terms  have  all 
been  shown  to  be  finite  after  integration  over  t'.  Similarly,  in 
the  terms  CT2  and  CT3,  the  T-factors  yield  convergent 
integrals  over  t'  (T3  behaves  just  like  T2),  hence  these 
terms  will  vanish  when  C  goes  to  zero. 
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The  only  terms  which  remain  are 
ALfct  =  -  C  cos (  6 ' —  6 )  T: 


and 

AL**  =  -  ?  Tj  . 

If  we  expand  cos ( 8 1  -  8),  the  second  term  is  proportional  to 
t2,  and  we  know  that  t2T],  is  integrable,  so  that  in  the 
limit,  we  need  retain  only 

ALfct  =  AL*'*’  =  -  ?  T1  . 


For  Ti  we  can  take  the  first  term  of  its  expansion; 


T1  = 


2*ax' 


since  all  others  are  finite  after  integration  over  t'.  The 
contribution  of  the  new  terms  thus  reduces  to 


Lfct( t) 


Kt(t' 


)  P' 


d  t 


Lim 


r 


Kt(t) 

2  IT 


Lim 

C+O 


p'dT 

2 


a  x 


with  an  exactly  analogous  expression  for  the  $-term.  In  the 
present  context,  we  may  think  of  6^  and  62  as  infinitesimal 
Using  the  expansions  derived  in  Section  2,  we  have 


£i  =  f  P*  \ 1/2 
a  P  +  C  cos 


1  + 


Tp 


(sin8  t 


cos  8  C )  +  . 


and 


2 

x  = 


(1 


2  2 

8  S)  tz  +  + 


where,  as  usual,  terms  of  higher  order  than  those  explictly 
indicated  turn  out  to  make  no  contribution  in  the  limit. 

We  define 


1  -  8  c 


and  substitute: 


Lim  C  j 
5*0 


P '  d  t 
2 

a  x 


Since 


T 


dt 


+  a 


is  finite  for  0=0,  the  corresponding  term,  when  multiplied  by 
C,  will  vanish  in  the  limit.  We  are  left  with 


Lim 

£-0 


1  -  B  c 


cos  S 
2  P 


O 


dt 


2  2 
T  +  O 


=  Lim 
o-*-0 


a 


2 


=  Lim 
o+O 


{ tan 


+  tan 


=  TT 


If  C  had  approached  the  surface  from  the  inside,  0  would  have 
passed  through  negative  values  to  zero,  and  the  result  would  have 
been  This  accounts  for  the  discontinuity  of  the  surface 

integral.  Altogether,  we  have  found 


t 


Lim  j  A  Ltt:(T)  Kt(t1)  P’  dT  =  -  i  Kfc(t) 
Z+o  -&l 


with  the  same  result  for  the  integral  of  AL^.  The  modal 
integral  equations  then  revert  to  the  forms  derived  from  the  Maue 
equation  directly. 
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